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Abstract 

Bose-Einstein correlations of pairs of identical charged pions produced in ha- 
dronic Z decays are analyzed in terms of various parametrizations. A good descrip- 
tion is achieved using a Levy stable distribution in conjunction with a model where 
a particle's momentum is correlated with its space-time point of production, the 
r-model. Using this description and the measured rapidity and transverse momen- 
tum distributions, the space-time evolution of particle emission in two-jet events is 
reconstructed. However, the elongation of the particle emission region previously 
observed is not accommodated in the r-model, and this is investigated using an ad 
hoc modification. 
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1 Introduction 



In particle and nuclear physics, intensity interferometry provides a direct experimental method 
for the determination of sizes, shapes and lifetimes of particle-emitting sources (for reviews 
see |TH5|) - In particular, boson interferometry provides a powerful tool for the investigation 
of the space-time structure of particle production processes, since Bose-Einstein correlations 
(BEC) of two identical bosons, first observed in 1959 [6]T7] and most recently at the Large Hadron 
Collider [5HT2]. reflect both geometrical and dynamical properties of the particle-radiating 
source. 

In e + e~ annihilation BEC have been observed [IB] to be maximal when the invariant mo- 
mentum difference of the bosons, Q = y/—(pi — P2) 2 , is small, even when one of the relative 
momentum components is large. This is not the case either in hadron-hadron interactions [H] 
or in heavy-ion interactions P2JH6], where BEC are found not to depend simply on Q, but 
to decrease when any of the relative momentum components is large, a behavior that can be 
described by hydrodynamical models of the source [5][T7]. 

The size (radius) of the source in heavy-ion collisions has been found to decrease with 
increasing transverse momentum, p t , or transverse mass, m t = \J 'm? + pf, of the bosons. This 
effect can also be explained by hydrodynamical models [17J[I8]. A similar effect has been seen 
in pp collisions [19], as well as in e + e~ annihilation [20Tf22] . 

A model for BEC in e + e~ annihilation, known as the r-model, which results in BEC de- 
pending simply on Q, rather than on its components separately, has been proposed [23J. In this 
model the simple Q dependence is a consequence of a strong correlation between a particle's 
four-momentum and its space-time point of production. This strong correlation also leads to a 
decrease of the observed source size with increasing m t |24j . 

On the other hand, studies of BEC in e + e~ annihilation at LEP have found that the BEC 
correlation function does depend on components of Q, the shape of the source being elongated 
along the event axis [2"2"ll2"5T - |2"8] . At HERA a similar elongation is observed in neutral current 
deep inelastic ep scattering [2H]. Note that here "source" refers not to the entire volume in 
which particles are emitted, but rather the smaller "region of homogeneity" from which pions 
are emitted that have momenta similar enough to interfere and contribute to the correlation 
function. The size of this region of homogeneity is sometimes referred to as the correlation 
length. 

The main purpose of the present paper is to test the r-model, but the question of to what 
extent BEC in e + e~ annihilation depend differently on different components of Q is also consid- 
ered. Here we study BEC in hadronic Z decay. We investigate various static parametrizations 
in terms of Q. Fitting over a larger Q range than in previous studies, we find that none of 
these parametrizations gives an adequate description of the Bose-Einstein correlation function. 
Next we investigate the r-model and find that it provides a good description of BEC, both 
for two-jet and three-jet events. The results for two-jet events, from fits in terms of Q and 
the transverse masses of the two pions with respect to the event axis, are used to reconstruct 
the complete space-time picture of the particle emitting source of two-jet events in hadronic Z 
decay. Finally, the relevance of the elongation mentioned above is investigated using an ad hoc 
modification of the r-model. 
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2 Analysis 



2.1 Data 

The data used in the analysis were collected by the L3 detector [3"0H3~4] at an e + e~ center-of- 
mass energy, y/s, of about 91.2 GeV. Calorimeter clusters having energy greater than 0.1 GeV 
are used to determine the event thrust [35] axis and to classify events as two- or three-jet events, 
which are analyzed separately, since differences in the space-time structure (and hence in the 
BEC parameters) could be expected, and indeed have been observed [26|l36]. 

Backgrounds to hadronic Z decays such as leptonic Z decays, beam-wall and beam-gas 
interactions, and two-photon events are excluded by requiring that the visible energy, E V i S , be 
within 0.5i/i and 1.5^/s, that the transverse and longitudinal energy imbalances be less than 
0.6E v[s and 0AE vis , respectively, and that the number of calorimeter clusters be at least 15. In 
order to ensure well-measured charged tracks, events are required to have the thrust direction 
within the barrel region of the calorimeters and the acceptance region of the tracking chamber: 
|cos(G)| < 0.74, where G is the angle between the thrust axis and the beam direction. High 
precision charged tracks are selected by requiring 

• transverse momentum greater than 150 MeV, 

• at least one hit in the inner region of the tracking chamber (TEC), 

• more than 25 hits in the entire TEC spanning at least 40 wires of the possible 62, 

• a distance, in the plane transverse to the beam, of closest approach to the interaction 
vertex less than 10 mm. 

Further, tracks lying in two small regions of azimuthal angle having less precise calibration 
are rejected: 45°-52.5° and 225°-232.5°. Since the resolution of the opening angle between 
pairs of tracks is crucial for the study of BEC, tracks are also rejected if there is no hit in the 
Z-chamber of the TEC. To further reject t + t~ events, the second largest angle, 02, between 
any two neighboring tracks in the transverse plane is required to be in the range 20°-170°. To 
further ensure that the event is well contained within the acceptance of the TEC, the thrust axis 
is determined using only the high precision tracks, and the event is rejected if |cos(Gtec)| > 0.7, 
where Gtec is the angle between this thrust axis and the beam. 

In total about 0.8 million events with an average number of about 12 high precision charged 
tracks are selected. This results in approximately 36 million like-sign pairs of charged tracks. 

The number of jets in an event is determined using the Durham jet algorithm [3"TH3"9"] with a 
jet resolution parameter y cut = 0.006, yielding about 0.5 million two-jet events and 0.3 million 
events having more than two jets. Since there are few events with more than three jets, the 
category of events with more than two jets is referred to as the three-jet sample. 

2.2 Bose-Einstein Correlation Function 

The two-particle correlation function of two particles with four-momenta p\ and pi is given by 
the ratio of the two-particle number density, p2{pi,P2), to the product of the two single-particle 
number densities, pi(pi) P\{p2) ■ Since we are here interested only in the correlation R2 due to 
Bose-Einstein interference, the product of single-particle densities is replaced by Po(pi,P2), the 
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two-particle density that would occur in the absence of Bose-Einstein correlations: 



fla(pi,P2) 



Po(Pl,P2) ' 



(1) 



This p 2 is corrected for detector acceptance and efficiency using Monte Carlo events generated 
by the Jetset Monte Carlo generator |3Q] with the so-called BE simulation of BEC [H] to 
which a full detector simulation has been applied [12], by multiplying the measured p 2 , on a 
bin-by-bin basis, by the ratio of p 2 of the generated events to p 2 of the generated events after 
detector simulation. For the detector-simulated events, as for the data, all charged tracks are 
used, whereas for the generator-level events only charged pions are used, since all measured 
tracks are assumed to be pions in the calculation of Q. Thus the Monte Carlo generator is used 
to extract the Q distribution for equally charged pion pairs from the observed Q distribution 
of all equally charged particle pairs. 

An event mixing technique is used to construct po, whereby all tracks of each data event 
are replaced by tracks from different events having a multiplicity similar to that of the original 
event. This is accomplished by first rotating each event to a frame whose axes are the thrust, 
major, and minor [33] directions and then assigning the events to classes based on the track 
multiplicity. Each multiplicity defines a class except for very low and very high multiplicities, 
in which case several adjacent multiplicities are grouped into the same class in order to obtain 
sufficient statistics. The replacement track is randomly chosen from tracks of the same sign in 
a randomly chosen event either in the same class as the data event or a neighboring or next- 
to-neighboring class. Ideally one would define the classes based on the total particle (charged 
plus neutral) multiplicity, since the aim is to use events where the available phase space of the 
particles is the same. However, the ratio of neutral to charged multiplicities is not constant 
and both multiplicities are subject to detector efficiency losses. In an attempt to take this into 
account we therefore also use events of nearby classes. 

A correction for detector acceptance and efficiency is applied to po m the same way as 
to p 2 . The mixing technique removes all correlations, e.g., resonances and energy- momentum 
conservation, not just Bose-Einstein correlations. Hence, po is also corrected for this by a 
multiplicative factor which is the ratio of the densities of events to mixed events found using 
events generated by Jetset without BEC simulation. 

Including all corrections, i? 2 is measured by 



where data, gen, det, gen-noBE refer, respectively, to the data sample, a generator-level Monte 
Carlo sample, the same Monte Carlo sample passed through detector simulation and subjected 
to the same selection procedure as the data, and a generator-level sample of a Monte Carlo 
generated without BEC simulation. 

We note that _R 2 is unaffected by single-particle acceptances and efficiencies. Being the same 
for p2 and po, these cancel in the ratio. Systematic uncertainties on R 2 are highly correlated 
point-to-point. Hence, in performing fits to i? 2 , only statistical uncertainties will be used in 
calculating the \ 2 which is minimized. 

In principle i? 2 should also be corrected for final-state interactions (FSI), both Coulomb 
and strong, or alternatively, the functional form used to fit i? 2 (Q) should be adapted to take 
account of FSI [331 135] . Coulomb interactions, being repulsive for like-sign charged pions, serve 
to increase the measured values of Q. This is often taken into account by weighting the pion 



R2 — (R2 data-R2 gen) / (-^2 dct-R2 gcn-noBE) , 



(2) 
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pairs with the inverse of the so-called Gamow factor [TJ. However, this factor strongly over- 
compensates, since it is derived assuming that the source is point-like and that all particles are 
pions [10]. Further the strong tt-tt FSI at short distances are attractive, further reducing the 
effect [17]. The net effect is small, amounting to about 3% at Q = [IS], comparable to the 
statistical uncertainty in the first bin of ^(Q) (c/. Figs.[2H3D. Further, if R 2 data is weighted 
to account for FSI, R2 gC n, RiAct and -R2gen-noBE need to be similarly weighted. This is because 
the Monte Carlo programs are tuned to data, and this tuning presumably accounts, at least 
partially, for the effects of FSI, even though the Monte Carlo program does not explicitly do 
so. The result is that the various weights, to a large degree, cancel. We have used unlike-sign 
charged pions to check that applying Gamow factors is not needed in the present analysis. 
The uncertainty on FSI is included in the systematic uncertaities by varying the fit range (cf. 
Section EQ)) . 



3 Parametrizations of BEC 
3.1 Dependence on Q 

Since the four-momenta of the two particles are on the mass-shell, the correlation function 
is defined in six-dimensional momentum space. Since BEC can be large only at small four- 
momentum difference Q, they are often parametrized in this one-dimensional distance measure. 
With a few assumptions [2HS1I7], the two-particle correlation function, Eq. (JT]), is related to the 
Fourier transformed source distribution: 



MQ) = 7 



l + A 



Kq) 



l + eQ), (3) 



where f(x) is the (space-time) density distribution of the source, and f(Q) is the Fourier 
transform of f(x). The parameter A is introduced as a measure of the strength of the correlation. 
It may be less than unity for a variety of reasons such as inclusion of misidentified non-identical 
particles or the presence of long-lived resonance decays if the particle emission consists of a 
small, resolvable core and a halo having experimentally unresolvable large length scales [Hi 
14*9] . The parameter 7 and the (1 + eQ) term parametrize possible long-range correlations not 
adequately accounted for in the reference sample. While there is no guarantee that (1 + eQ) 
is the correct form, we will see that it does provide a good description of R2 in the region 
Q > 1.5 GeV. In fact, if the reference sample indeed removes only BEC, we expect R2 to 
contain no long-range correlations and hence to find e = 0. 

There is no reason, however, to expect the hadron source to be spherically symmetric in jet 
fragmentation. Recent investigations have, in fact, found an elongation of the source along the 
jet axis [221E5H2E] • While this effect is well established, the elongation is actually only about 
20%, which suggests that a parametrization in terms of the single variable Q, may be a good 
approximation. 

This is not the case in heavy-ion and hadron-hadron interactions, where BEC are found 
not to depend simply on Q, but on components of the four-momentum difference separately 
[5*l [T4TfT7t [r5] . However, in e + e~ annihilation at lower energy [13] it has been observed that Q is 
the appropriate variable. We checked this [SD] both for all and for two-jet events: We observe 
that i? 2 does not decrease when both q 2 = (pi — p 2 ) 2 and q$ = {E\ — E 2 ) 2 are large while 
Q 2 = q 2 — q 2 is small, but is maximal for Q 2 = q 2 — q 2 = 0, independent of the individual 
values of q and q . In a different decomposition, Q 2 = Q\ + Ql E , where Q\ = (p t i — A2) 2 is 
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Figure 1: (a) R 2 for two-jet events as function of the squares of the transverse momentum 
difference, Q 2 = (p t i — Pt2) 2 -, and the combination of longitudinal momentum difference and 



energy difference, Ql E = (p u - pi 2 ) - (E 1 - E 2 ) . (b) R 2 vs. Q\ when Q\ 
corresponds to Q 2 pa 0. 
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the component transverse to the thrust axis and Ql E = (Pu — P12) 2 — {E\ — E 2 ) 2 combines the 
longitudinal momentum and energy differences, we find R 2 to be maximal along the line Q = 0, 
as is shown in Fig.[U and fits, though of poor x 2 (confidence level of about 1%), are consistent 
with equal radii, as has previously been observed by aleph [5Tj . 

We conclude that a parametrization in terms of Q can be considered a reasonable ap- 
proximation for the purposes of this article. We shall return to the question of elongation in 
Section 



3.2 Symmetric parametrizations 

The simplest assumption is that the source has a symmetric Gaussian distribution with mean 
fi = and standard deviation R. In this case f(Q) = exp (i^iQ — ^ RC ^ \ and 

RziQ) = 7 [l + A exp (-(RQ) 2 )] (1 + eQ) . (4) 

However, this parametrization is often found to be inadequate in its description of data. 

A model-independent way to study deviations from the Gaussian parametrization is to 
use J3JE5JE3] the Edgeworth expansion [51] about a Gaussian. Keeping only the first non- 
Gaussian term, we have 

R*{Q) = 7 (l + A exp {-(RQ) 2 ) [l + ^ 3 (JJQ)]) (1 + eQ) , (5) 

where k is the third-order cumulant moment and H^(RQ) = (\/2RQ) 3 — 3\/2RQ is the third- 
order Hermite polynomial. Note that the second-order cumulant corresponds to the radius R. 
Eq. d5j reduces to Eq. fflh for k = 0. 
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Another way to depart from the assumption of a Gaussian source is to replace the Gaussian 
by a symmetric Levy stable distribution, which is characterized by three parameters: xq, R, 
and a. Its Fourier transform, f(Q), has the following form: 



f(Q) = exp ( iQx 



\RQ\ 



2 



a 



) 



(6) 



where the index of stability, a, satisfies the inequality < a < 2. The case a = 2 corresponds 
to a Gaussian source distribution with mean xq and standard deviation R. For more details, 
see, e.g., [55]. Then R2 has the following form [56] : 




These three parametrizations are fitted to the data for both two- and three-jet events. Fits 
of the Gaussian parametrization, Eq. (j4]), to the data (shown in Figs. [2k, and [3^ for two- and 
three-jet events, respectively) result in unacceptably low confidence levels: 10~ 15 (x 2 = 246 for 
96 degrees of freedom) for two-jet and much worse (x 2 = 456) for three-jet events. The fit is 
particularly bad at low values of Q, where both for two- and three-jet events R2 is much steeper 
than a Gaussian. 

A fit of the Edgeworth parametrization, Eq. (JS]), to the two-jet data, shown in Fig.[2b, finds 
k = 0.74 ± 0.07, about 10 standard deviations from the Gaussian value of zero. The confidence 
level is indeed much better than that of the purely Gaussian fit, but is still poor, approximately 
10 -5 . Close inspection of the figure shows that the fit curve is systematically above the data 
in the region 0.6-1.2 GeV and that the data for Q > 1.5 GeV appear flatter than the curve, as 
is also the case for the purely Gaussian fit. Similar behavior, but with a worse x 2 ; is observed 
for three-jet events (Fig. [3b). 

The fit of the Levy parametrization, Eq. ([?]), to the two-jet data, shown in Fig. [2b, finds 
a = 1.44 ± 0.06, far from the Gaussian value of 2. The confidence level, although improved 
compared to the fit of Eq. (J4[) is still unacceptably low, approximately 10~ 8 . Similar behavior, 
with a worse y 2 , is seen for three-jet events (Fig.[3b). 

Both the symmetric Levy parametrization and the Edgeworth parametrization do a fair job 
of describing the region Q < 0.6 GeV, but fail at higher Q. R2 in the region Q > 1.5 GeV 
is nearly constant (~ 1). However, in the region 0.6-1.5 GeV R2 has a smaller value, dipping 
below unit}0, which is indicative of an anti-correlation. This is clearly seen in Figs. [2] and [3] 
by comparing the data in this region to an extrapolation of a linear fit, Eq. ([4]) with A = 0, in 
the region Q > 1.5 GeV. The inability to describe this dip in R2 is the primary reason for the 
failure of both the Edgeworth and symmetric Levy parametrizations. Failure to describe this 
anti-correlation region, i.e., the dip below unity, is an inherent feature of any parametrization 
of the form Eq. ([3]). We note that this dip is less apparent if one only plots (and fits) R2 for 
Q < 2 GeV as has usually been done in the past, the dip being accommodated by a relatively 
large value of e. 

3.3 Time dependence of the source 

The parametrizations discussed so far, which have proved insufficient to describe the BEC, 
all assume a static source. The parameter R is a constant. It has, however, been observed 
that R depends on the transverse mass, m t = m 2 + p 2 = sjE 2 — p 2 , of the pions [20H22] . 

*More correctly, dipping below the value of 7(1 + e). 





R 2 (Q) = 7 [1 + A exp (-0RQ) Q )] (1 + eQ) . 



(7) 
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Figure 2: The Bose-Einstein correlation function R 2 for two-jet events with the results of fits of (a) the Gaussian, Eq. 01]), (b) the 
Edgeworth, Eq. (J3J), and (c) the symmetric Levy, Eq. (J7J), parametrizations. Also plotted is A, the difference between the fit and 
the data. The dashed line represents the long-range part of the fit, i.e., 7(1 + eQ). The dotted line represents a fit with A = to 
Q > 1.52 GeV. 
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Figure 3: The Bose-Einstein correlation function R 2 for three-jet events with the results of fits of (a) the Gaussian, Eq. fl4]), (b) the 
Edgeworth, Eq. ([5]) , and (c) the symmetric Levy, Eq. (J7J), parametrizations. Also plotted is A, the difference between the fit and 
the data. The dashed line represents the long-range part of the fit, i.e., 7(1 + eQ). The dotted line represents a fit with A = to 
Q > 1.52 GeV. 



It has been shown [53 [58] that this dependence can be understood if the produced pions 
satisfy, approximately, the (generalized) Bjorken-Gottfried condition J59HM], whereby the four- 
momentum of a produced particle and the space-time position at which it is produced are 
linearly related. Such a correlation between space-time and momentum-energy is also a feature 
of the Lund string model, which, as incorporated in Jetset, is very successful in describing 
detailed features of the hadronic final states of e + e~ annihilation. Not only inclusive single- 
particle distributions and event shapes, but also correlations [65H67] are well described. 



3.3.1 The r-model 

In section 13.11 we have seen that BEC depend, at least approximately, only on Q and not on 
its components separately. This is a non-trivial result. For a hydrodynamical type of source, 
on the contrary, BEC decrease when any of the relative momentum components is large OUT] . 
Further, we have seen that R2 shows anti-correlations in the region 0.6-1.5 GeV, dipping below 
its values at higher Q. 

A model which predicts such Q-dependence, as well as the absence of dependence on the 
components of Q separately, is the so-called r-model (23]. Further it incorporates the Bjorken- 
Gottfried condition and predicts a specific transverse mass dependence of R2, which we subject 
to an experimental test here. 

In this model, it is assumed that the average production point in the overall center-of-mass 
system, x = (t, r x , r y ,f z ), of particles with a given four-momentum p = (E,p x ,p y ,p z ) is given 
by 

z"(p") = arp^ . (8) 



In the case of two-jet events, a = l/m t where m t is the transverse mass and r = yt — r z 2 is 
the longitudinal proper timeo For isotropically distributed particle production, the transverse 

mass is replaced by the mass in the definition of a and r by the proper time, \Jt 2 — r^ — r y — r z 2 . 
In the case of three-jet events the relation is more complicated. 

The second assumption is that the distribution of x M (p M ) about its average, S^x^ip^) — 
x^ijf)), is narrower than the proper-time distribution, H(t). Then the two-particle Bose- 
Einstein correlation function is indeed found to depend on the invariant relative momentum Q, 
rather than on its separate components, as well as on the values of a of the two particles [23] : 

i? 2 (p^ 2 ) = l + Re^(^)if(^) , (9) 

where H(u) = J drH(r) exp(iur) is the Fourier transform (characteristic function) of H{t). 
Note that H(t) is normalized to unity. 

Since there is no particle production before the onset of the collision, H{j) should be 
a one-sided distribution. In the leading log approximation of QCD the parton shower is a 
fractal [58TT70] . Further, a Levy distribution arises naturally from a fractal [71]. We are thus 
led to choose a one-sided Levy distribution for H(t) [23] • The characteristic function of H(t) 



*The terminology 'longitudinal' proper time and 'transverse' mass seems customary in the literature even though 
their definitions arc analogous t = \Jt 2 — r z 2 and m t = \J E 2 — p% . 
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can then be written [56] (for q^I) as3 

1 



H u) 



exp 



(Ar|w|) Q (l - zsign(w)tan ) 



(10) 



where the parameter To is the proper time of the onset of particle production and At is a 
measure of the width of the proper-time distribution. Using this characteristic function in 
Eq. 09]), and incorporating the factor A and the long-range parametrization, yields 



R 2 (Q, ai, a 2 ) = 7 \ 1 + Acos 



T Q 2 (a 1 + a 2 ) 



■ exp 



ArQ 2 \ a a« + a% 



(11) 



Note that the cosine factor generates oscillations corresponding to alternating correlated 
and anti-correlated regions, a feature clearly seen in the data (Figs. [2] and [3]). Note also that 
since a = l/m t for two-jet events, the r-model predicts a decreasing effective source size with 
increasing m t . These features are subjected to a quantitative test in the following sections. 



3.3.2 The r-model for average a 

Fits of Eq. ffTT]) are difficult since i? 2 depends on three variables: Q, a\ and a 2 . Further, we have 
a simple expression for a only for two-jet events. Therefore, we first consider a simplification of 
Eq. (TTT]) obtained by assuming (a) that particle production starts immediately, i.e., r = 0, and 
(b) an average a-dependence, which is implemented by introducing an effective radius defined 
by 

'Axro? + a| (12) 



R 



2a 



This results in 



R 2 (Q) = 7 [1 + Acos ((R*Q) 2a ) exp (- (RQ) 2a )] (1 + eQ) 



where R a is related to R by 



R 



2a 



fair 
tan — 
V 2 



R 



2a 



(13) 



(14) 



Fits of Eq. ( JT3"j) are first performed with R a as a free parameter. The fit results obtained for 
two- and three-jet events are listed in Table [T] and shown in Fig.H] for two-jet events and in 
Fig. [5] for three-jet events. They have acceptable confidence levels, describing well the dip below 
unity in the 0.6-1.5 GeV region, as well as the peak at low values of Q. 

As shown in Table [T] the estimates of some fit parameters are rather highly correlated. 
Taking these correlations into account, the fit parameters for the two-jet events satisfy Eq. f)14p . 
the difference between the left- and right-hand sides of the equation being less than a standard 
deviation. For three-jet events the difference amounts to about 1.5 standard deviations. 

Note that no significant long-range correlation is observed: e is zero within 1 standard 
deviation, and fits with e fixed to zero find the same values, within 1 standard deviation, of the 
parameters as the fits of Table [T] Thus the method to remove non-Bose-Einstein correlations 
in the reference sample is apparently successful. 



*Our notation, H(uj), corresponds to Nolan's [5S] S(a, (3, 7, 8; 1) parametrization with (3 
and 5 tq . The special case of a = 1 is also given by Nolan [SJj . 
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Figure 4: The Bose-Einstein correlation function R2 for two-jet events. The curve corresponds to the fit of the one-sided Levy 
parametrization, Eq. ( I13j) . with the parameter R. A (a) free and (b) constrained by Eq. (I14j) . The results of the fits are given in Tables [U 
and [21 respectively. Also plotted is A, the difference between the fit and the data. The dashed line represents the long-range part of 
the fit, i.e., 7(1 + eQ). 
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Figure 5: The Bose-Einstein correlation function R 2 for three-jet events. The curve corresponds to the fit of the one-sided Levy 
parametrization, Eq. (1131) . with the parameter R a (a) free and (b) constrained by Eq. (1141) . The results of the fits are given in Tables Q] 
and El respectively. Also plotted is A, the difference between the fit and the data. The dashed line represents the long-range part of 
the fit, i.e., 7(1 + eQ). 



Table 1: Results of fits of Eq. ( f!3|) for two-jet and three-jet events, which are shown in Figs.Hk 
and|5k, respectively. The first uncertainty is statistical, the second systematic (see Section f3T4|) . 



parameter two-jet three-jet 

A 

a 

R (fm) 
i?a (fm) 
e (GeV" 1 ) 

7 

X 2 /DoF 91/94 84/94 

confidence level 57% 76% 



correlation coefficients two-jet three-jet 



a 


R 




e 


7 


a 


R 




e 


7 


A -0.72 


0.95 


0.52 


-0.22 


0.27 


A -0.78 


0.97 


0.65 


-0.32 


0.38 


a 


-0.62 


-0.92 


0.72 


-0.79 


a 


-0.77 


-0.95 


0.74 


-0.81 


R 




0.38 


-0.10 


0.13 


R 




0.61 


-0.27 


0.32 








-0.89 


0.94 








-0.88 


0.93 










-0.97 


C 








-0.98 



0.63 ± 0.031^ 
0.41 ± 0.02+°;°^ 
0.79 ± 0.04±g;?| 



0.69 ± 0.041°;^ 
0.001 ± 0.002; 



-0.005 
-0.008 



0.988 ± 0.005lg;g^ 



0.92 ± 0.051^ 
0.35 ± 0.01+°;^ 
1.06 ± 0.05^|? 
0.85 ± 0.04ig;^ 
0.000 ± 0.00218X 
0.997 ± 0.005±S;gi| 



Fit results imposing Eq. (I14p are given in Table [2j For two-jet events, the values of the 
parameters are comparable to those with i? a free. For three-jet events, the imposition of 
Eq. fll4p results in values of a and R closer to those for two-jet events. Unlike the fits with R a 
free, the values of e differ somewhat from zero, which could indicate a slight deficiency in the 
description of BEC in the same way, though on a much smaller scale, as seen in the Edgeworth 
and symmetric Levy fits of Figs. [2] and [3j 

3.3.3 The r-model with ra t dependence 

For two-jet events, a = l/m t , while for three-jet events the situation is more complicated. We 
therefore limit fits of Eq. ffTTT) to the two-jet data. For each bin in Q the average values of 
m tl and m t2 are calculated, where m tl and m t2 are the transverse masses of the two particles 
making up a pair, requiring m tl > m t2 . Using these averages, Eq. ffTT]) is fit to i? 2 (Q)- The 
fit results in To = 0.00 ± 0.02 fm, and we re-fit with r fixed to zero. The results are shown 
in Fig.[6]and Table [31 The parameters a, At and A are highly correlated. Eq. (TlT]) for two-jet 
events can be simplified by requiring m t i ~ m t2 . This leads to fit results [50] consistent with 
those without this simplification. 

Since the r-model describes the m t dependence of i? 2 , an important test of the model is 
that its parameters, a, At, and r , not depend on m t . Note that the parameter A, which is not 
a parameter of the r-model, can depend on m t , e.g., as a result of resonances |48,49j. The large 
correlations between the fit estimates of A and those of the r-model parameters complicate the 
testing of m t -independence. We perform fits in various regions of the m tl -m t2 plane keeping a 
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Table 2: Results of fits of Eq. (1131) imposing Eq. (j!4p for two-jet and three-jet events, which are 
shown in Figs. lib and[5b, respectively. The first uncertainty is statistical, the second systematic 
(see Section 13. 4p . 





parameter 


two-jet 






three-jet 
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0.61 ± 0.03+°;°® 




0.8- 


i ± 0.04±g;g* 
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0.44 ± 0.01±g;g| 




0.42 ± 0.01±g;gl 




R (fm) 


0.78 ± O.OAtofe 




0.9; 


3 ± 0.04+°;^ 




e (GeV^ 1 ) 


0.005 ±0.001 ±0.003 


o.oo, 


3 ± 0.001 ±0.005 
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0.979 ± 0.002tg:oos 


0.977 ± O.OOli^s 




x 2 /DoF 


95/95 






113/95 




confidence level 


49% 






10% 


correlation coefficients 


two- jet 






three-jet 




a R 


e 7 




a 


R e 7 


A 


-0.92 0.96 


0.29 -0.34 


A 


-0.93 


0.96 0.24 -0.27 


a 


-0.96 


-0.35 0.42 


a 




-0.97 -0.31 0.37 


R 




0.21 -0.26 


R 




0.20 -0.25 


e 




-0.89 


c 




-0.89 



Table 3: Results of the fit of Eq. (II ip for two-jet events, which is shown in Fig.0 The parameter 
T is fixed to zero. The first uncertainty is statistical, the second systematic (see Section E3J • 



parameter 



A 0.58 ± 0.03tg;2f 

a 0.47 ± 0.01+°,;°^ 

At (fm) 1.56 ±0.12^;| 

e (GeV" 1 ) 0.001 ± 0.001 ±0.003 

7 0.988 ± 0.002+°;°°^ 
X 2 /DoF 90/95 
confidence level 62% 
correlation coefficients 
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At 


e 
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A -0.92 


0.95 


0.30 


-0.35 


a 


-0.96 


-0.37 


0.45 


At 




0.23 


-0.29 


c 






-0.89 



11 



1.8 



1.6 



1.4 - 



CV2 



1.2 



1.0 - 



0.8 
0.02 
< 0.00 
-0.02 
-0.04 



L3 

two— jet 
full t model 








Q (GeV) 



Figure 6: The Bose-Einstein correlation function R 2 for two-jet events. The curve corresponds 
to the fit of Eq. (Hip . The results of this fit are given in Table Also plotted is A, the difference 
between the fit and the data. The dashed line represents the long-range part of the fit, i.e., 
7(1 + eQ). 
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Table 4: Confidence levels and the values of A found in fits of Eq. (11 ip for two-jet events in 
various regions of the m t i-m t 2 plane with a and At fixed to the result of the fit to the entire 
plane. 



m t regions (GeV) 


averaj 




POD fl Pi PT) PP 
V^VJllllVJ-CllV^C 




77? 4-1 
"HI 




m t (GeV) 
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1C V C± 
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n 1 q 


n 1 q 
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n on _|_ n no 
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— n 
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0.37 


74 


n QO _)_ n no 
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0.52 


1 3 


n Q7 + n (14 


0.26 - 0.42 


0.14 
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0.25 


0.26 


22 


0.53 ± 0.02 
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0.32 


0.33 


33 


0.80 ±0.03 


0.46 - 0.58 


0.30 


- 0.42 


0.43 


0.44 


34 


0.91 ±0.04 
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0.42 


-4.14 


0.65 


0.65 


66 


1.01 ±0.05 


0.42 - 0.62 


0.14 


- 0.22 


0.34 


0.34 


17 


0.41 ±0.03 


0.46 - 0.70 


0.22 


- 0.30 


0.41 


0.41 


55 


0.64 ±0.03 


0.58 - 0.82 


0.30 


- 0.42 


0.52 


0.52 


59 


0.70 ±0.04 


0.86 - 1.22 


0.42 


- 4.14 


0.80 


0.81 


24 


0.66 ±0.05 


0.70 - 4.14 


0.22 


- 0.30 


0.59 


0.65 


4 


0.37 ±0.04 


0.82 - 4.14 


0.30 


- 0.42 


0.71 


0.76 


11 


0.56 ±0.05 



and At fixed at the values obtained in the fit to the entire m t plane (Table [3]). The regions 
are chosen such that the numbers of pairs of particles in the regions are comparable. The m t 
regions and the confidence levels of the fits are shown in Table S] along with the values of A. The 
confidence levels are seen to be quite reasonably distributed in agreement with the hypothesis 
of m t -independence of the parameters of the r-model (a, At, and To). 

3.4 Systematic Uncertainties 

The following sources of systematic uncertainty are investigated: 

• Event and track selection: The cuts used in selecting hadronic events and high precision 
tracks are varied within reasonable limits: (a) E vis within yfs ± 0.45-y/i to within yfs ± 
0.55-^/i; (b) transverse energy imbalance from 0.5-E v i S to 0.7-E v j S and longitudinal energy 
imbalance from 0.3-E v i S to 0.5-E v i S ; (c) minimum number of calorimeter clusters from 13 to 
17; (d) maximum value of |cos(0)| from 0.707 to 0.777 and maximum value of |cos(6tec)| 
from 0.662 to 0.736; (e) transverse momentum of a track from 100 MeV to 200 MeV; (f) 
minimum of 21 hits spanning at least 48 wires to 31 hits spanning at least 32 wires; (g) 
maximum distance of closest approach from 5 mm to 15 mm. 

• Fit range: The fits presented above were performed in the range < Q < 4 GeV. The 
fits were repeated in the ranges 0.04 < Q < 4 GeV, < Q < 3 GeV and < Q < 5 GeV. 
Note that uncertainty on the contribution of FSI is included by considering the range 
0.04 < Q < 4 GeV, FSI being most important at the smallest values of Q. 
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• Mixing: The event mixing method to construct po uses tracks from events having similar 
multiplicity. The definition of "similar" was varied by changing the range of multiplicity 
classes from which the mixed tracks are chosen from to 8, i.e., from demanding the 
same multiplicity class as the data event to including classes within ±4 of that of the 
data event. 

• Monte Carlo: Events generated by Jetset with no BEC simulation or by Herwig [72] 
are used instead of Jetset with BEC for the correction of p2 for detector acceptance, 
efficiency and non-pion pair background and PYTHIA [73] and HERWIG [71] are useqj 
instead of Jetset without BEC to correct po- 

For each source the root mean square of the deviations from the values measured in the standard 
selection is taken as the systematic uncertainty, positive and negative deviations being treated 
separately. The systematic uncertainties from the four sources are added in quadrature to 
obtain the total systematic uncertainty. For the physically interesting parameters (a, R, R a , 
At), the largest contributions come from Monte Carlo variation and/or mixing. The other 
parameters sometimes receive significant contributions from fit range and/or track and event 
selection. 



4 The emission function of two-jet events 

Within the framework of the r- model, we now reconstruct the space-time picture of the emitting 
process for two-jet events. The emission function in space-time, S rr (f,T), normalized to unity, 
is given in the r-model by 



1 d 4 n 
S r Jr, r = — 



n drdr 



where n and n are, respectively, the number and average number of pions produced, and where 
p p (p) is the single-particle momentum distribution, 

= dJ^dS ■ (16) 

which is normalized to the mean multiplicity, n. Note that, given the r-model correlation 
between space-time and momentum space, 

m t = t = . (17) 

T 2 _ ( r 2 _|_ r 2) 

Using t 2 = t 2 — r 2 and Eq. (|T7I) , we can rewrite Eq. (I15H in terms of r and t: 

SA^t) = -^=-J It ^t)H(r=y^^)p p lp= , (18) 

ndtdr n J \ / t 2 _ (r 2 + -2 + r 2)) J 



* Herwig was modified to use the decay and BEC routines of Pythia. Fragmentation parameters of both 
Herwig and Pythia were tuned using L3 data with BEC simulated by the BE32 Gaussian algorithm [75] , 
The events are generated using these parameters but without BEC simulation. 
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Figure 7: The proper-time distribution, H(t), for a = 0.47, tq = and At = 1.56 fm. H(t) 
was calculated using the program STABLE |76j . 



is the Jacobian of the variable transformation. 

Given the symmetry of two-jet events, it is convenient to write 5* in cylindrical coordinates 
and average over the azimuthal angle. To simplify the reconstruction of S we assume that the 
momentum distribution p p can be factorized in the product of transverse and longitudinal dis- 
tributions, e.g., p yPt = Py(y)p Pt (pt), where y is the rapidity and p t is the transverse momentum. 
This is found to be a reasonable approximation except at very high values of \y\ or p t [50J. 

Using H{t) as obtained from the fit of Eq. (11 11) (Table [3]), which is shown in Fig. [71 together 
with the inclusive p t and rapidity distributions^ which are shown in Fig. [HI the full emission 
function is reconstructed. 

Integrating Eqs. (1151) and ffTHl) over the transverse coordinates results, respectively, in 



* These distributions of high precision charged tracks from L3 two-jet events are corrected for detector accep- 
tance and efficiency bin-by-bin by the ratio of the distributions of generator-level to detector-level Monte Carlo 
samples. For the purpose of the reconstruction, the pt and rapidity distributions have been parametrized 



as Ppt ( Pt ) = exp (-8.340p t + 1.249^ - 0.1050p?) (l9.52p t + 182. bp 2 t - lUApl + 90.37rf + 48.84p 5 t ) and 



where 




(19) 




(20) 



p y {y) = cxp ( - [(-0.6 
-0.2843 |y 5 | +0.07887?; 6 ) . 
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Figure 8: The (a) p t and (b) rapidity distributions of two-jet events and the parametrizations 
used in the reconstruction of the source function. 




Figure 9: The temporal-longitudinal part of the source function of two-jet events, (a) S(t),t), 
Eq. flU and (b) S(r z ,t), Eq. (J2TJ , normalized to unity. 

and 

**<■■■■ 4) 4^ - ^7P=j ff ( T = ^) * (» = " =5 b £S • (21) 

These are plotted in Fig.|9j They exhibit a "boomerang" shape with maxima at low values of 
r and rj or t and r z , but with tails reaching out to very large values, a feature also observed 
in hadron-proton [77] and heavy ion collisions [75]. Note however, that in the case of two-jet 
events, the emission function S^ir), r) has two maxima, as does S zt (r z ,t), as is seen in Fig.0 
These two different maxima correspond to the two-jet structure of the events. This is in marked 
contrast with the reconstructed emission function in hadron-proton reactions at ^/s = 22 GeV, 
where only a single maximum is observed [77] . 

The transverse part of the emission function is obtained by integrating Eq. ffl5|) over r z and 
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averaging over the azimuthal angle: 




Fig. ITOl shows the transverse part of the emission function, Eq. (122ft . for various proper times. 
Particle production starts immediately, increases rapidly and decreases slowly A ring-like 
structure, similar to the expanding, ring-like wave created by a pebble in a pond, is observed. 
These pictures together form a movie representing this temporal process, which takes place, to 
our knowledge, at the shortest time scale ever measured^ In the movie two ring-like structures 
separate along the thrust axis, their radii increasing. They thus span the surface of two back- 
to-back cones whose tips meet at the origin. The rings are very faint at the beginning, then 
quickly brighten as the particle production probability increases to a sharp peak at r m 0.035 fm 
(« 0.12 x 10~ 24 s). These separating and expanding rings then start to fade gradually with a 
characteristic duration of At 1.56 fm (« 5.2 x 10~ 24 s) following a power-law tail. Note 
however that the shape of these rings is in marked contrast to the more smoothly edged rings 
of fire found in hadron-proton reactions [77]. In our case, these rings are reconstructed from 
the Bose-Einstein correlation functions and from the measured p t and rapidity spectra of two- 
jet events without any reference to temperature or flow, indicating the non-thermal nature of 
particle production in e + e~ annihilation. 



5 Test of dependence of BEC on components of Q 

In this section we return to the question of a possible elongation of the region of homogeneity 
and the possible dependence of the Bose-Einstein correlation function on components of Q. 

The r-model predicts that the two-particle BEC correlation function R 2 depends on the 
two-particle momentum difference only through Q, i.e., that a single radius parameter applies 
to all components of Q. However, previous studies [2"2" ||2"5"H2"%] have found that different radii 
are required for the different components of Q in order to fit the data, and indeed that the 
shape of the region of homogeneity is elongated along the event (thrust) axis. The question is 
whether this is an artifact of the Edgeworth or Gaussian parametrizations used in these studies 
or shows a limitation of the r-model. 

5.1 Test of elongation and the r-model 

In the previous studies of elongation, Q 2 is split into three components, 

Q 2 = Ql + Qside + QL - ( A£) 2 (24) 

= Ql + Qldc + QL (1 - n , where = ^1±|^ (25) 

Mil + hi2 

*Animated gif files covering the first 0.15fm (0.5 x 10~ 24 s) are available [75] . 
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Figure 10: The transverse source function S(x,y,r), Eq. (122]) . normalized to unity, and its 
transverse profile for various proper times ranging from 0.01 fm to 0.15 fm. 
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in the Longitudinal Center of Mass System (LCMS) of the pair. The LCMS frame is defined 
as the frame, obtained by a Lorentz boost along the event axis, where the sum of the three- 
momenta of the two pions (p\ + P2) is perpendicular to the event axis. Assuming azimuthal 
symmetry about the event axis suggests that the region of homogeneity have an ellipsoidal 
shape with one axis, the longitudinal axis, along the event axis. The longitudinal axis would 
then be longer, shorter, or equal to the other two (transverse) axes, which are of equal length. 

In the LCMS frame the event axis is referred to as the longitudinal direction; the direction 
of p\ + P2 as the out direction; and the direction perpendicular to these as the side direction. 
The components of \pi — p 2 \ along these directions are denoted by Ql, Qout, and Q S id e - In the 
Gaussian and Edgeworth parametrizations R 2 Q 2 is then replaced in the fitting function by 

R 2 Q 2 =^ RlQl + -RsidcQsidc + PoutQout • (26) 



The longitudinal and transverse sizes of the source are measured by i?L and -R S id e , respectively, 
whereas the value of p out reflects both the transverse and temporal sizesQ Note that the LCMS 
and the rest frame of the pair differ only by a Lorentz boost of velocity along the out direction. 
Hence -Rl and -R S id c also measure the longitudinal and transverse size in the pair rest frame, 
while p out may be interpreted as an average of the transverse size in the rest frame boosted to 
the LCMS: Po 2 ut = r 2 ut (l - (5 2 ). 

In our previous analysis [25] we found, using all events, R side /R L = 0.80 ± 0.021°;°! and 
Rside/ Rl = 0.81 ± 0. 02l ) ; J for the Gaussian and Edgeworth parametrizations, respectively. 
Since the present analysis uses a somewhat different event sample and mixing algorithm, we 
have repeated these fits using the same binning of 0.08 GeV for 0.00 < Qi < 1.04 GeV (z=L, 
side, out). We find consistent values [50] of R S id e /Rh'- 0.76 for both parametrizations. We also 
find the elongation to be greater for two-jet events than for three-jet events, as has previously 
been observed by opal [26] • Using the Edgeworth parametrization the elongation is (statistical 
uncertainties only) R S ia e / Rl = 0.65 ± 0.02 for two-jet events, 0.76 ± 0.02 for all events, and 
0.84 ± 0.02 for three-jet events. 

We next investigate whether the r-model parametrizations could accommodate an elonga- 
tion. This is done by incorporating Eq. (I26T) in Eq. (TT5I) . If Eq. (FT4T) is imposed, this results 
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R2{Q)=i 



1 + A cos (tan (^) A 2a ^j exp (-A 2a ) 



1 + e^Qh + CsideQside + CoutQout) , (27) 



where 



A 2 



RlQl + -^sidoQsidc + PoutQout > 



(28) 



and where we use the same long-range parametrization as in the above Edgeworth and Gaussian 
parametrizations. We have attempted the long-range parametrization (1 + eQ), but that results 
in unacceptable x 2 . To extend the range of Q included in the fit, the bin size is increased to 
0.16 GeV for Q { > 0.88 GeV. Results for two-jet events of a fit of Eq. ([27]) for Q < 4 GeV are 
shown in Table [5j A clear preference for elongation (i? si d c /-RL < 1) is seen. The value of R 
found in the corresponding fit of i?2(Q) (Table [2]) lies between the values of -R S id c and R^ found 
here. Furthermore, the value of the elongation agrees well with that found above for two-jet 
events using the Edgeworth parametrization. 



*In the literature the coefficient of Q^, ut in Eq. is usually denoted R^, ut . We prefer to use p^ ut to emphasize 
that, unlike i?L and i? s ido, Pout contains a dependence on /3 and to differentiate it from i? ou t in Eq. (|31a|) below. 
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Table 5: Results for two-jet events of fits of Eq. ( 12 7p . Uncertainties are statistical only. 
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A 


n aq -+- n D9 


f\: 
IX 


n 4« _i_ n m 
u.^iu ziz u.ui 


Rt (fm) 
- n 'L \ L11L ) 


U.OJ _I_ U.Ut: 


-Kside/ -Kl 


n ^ 1 _i_ n no 
U.Ol ± U.U2 


Pout/ Rl 


0.66 ±0.02 


e L (GeV- 1 ) 


0.001 ±0.001 


e side (GeV- 1 ) 


-0.076 ± 0.003 


e ou t (GeV- 1 ) 


-0.029 ±0.002 


7 


1.011 ±0.002 


X 2 /DoF 


14847/14921 


CL 


66% 



5.2 Direct test of R 2 dependence on components of Q 

To directly test the hypothesis of no separate dependence on components of Q, we investigate 
two decompositions of Q: 

Q 2 - Qle ± Qsidc + Qlut (29a) 

Q 2 = Ql + QLc + dt (29b) 

where Ql E = Ql-(AE) 2 and q% nt = Ql nt — (AE) 2 . The decomposition of Eq. (I29ajl corresponds 
to the LCMS frame where the longitudinal and energy terms are combined; its three components 
of Q are invariant with respect to Lorentz boosts along the thrust axis. The decomposition 
of Eq. (129bp corresponds to the LCMS frame boosted to the rest frame of the pair; its three 
components of Q are invariant with respect to Lorentz boosts along the out direction. 

The test then consists of replacing R 2 Q 2 by R 2 times one of the above equations and 
comparing a fit where the coefficients of all three terms are constrained to be equal to a fit 
where each of these coefficients is a free parameter. 

Modifying Eq. (fT3"j) in this way, and imposing Eq. (I14p . results in 



R2(Q) = i 

where 



1 ± A cos ( tan ( 



(air 



B 2a ^j exp (-B 2a ) 



(30) 



B 2 — -RleQlE ± -^sideQside ± -^outQout (31a) 

6 = 1 + CleQlE + ^sideQside + CoutQout (31b) 

or 

B 2 = R\Q\ + -RgideQsidc + r out<?out (32a) 

b = 1 + ClQl + e sideQside + Cout?out • (32b) 

Note that, since AE vanishes in the rest frame of the pair, the quantity r out is a direct 

measure of the spatial extent of the source in this frame whereas p out in Eq. (1281) and i?LE in 
Eq. (13 laft are also sensitive to the temporal distribution of the source. 
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The results of the four fits of Eq. (1301) are shown in Table [6] for Q < 4 GeV. For both 
parametrizations the fit allowing separate dependence on the components of Q has a higher 
confidence level than the fit which does not. The fit using Eq. (1311) attains a confidence level of 
2% when no separate dependence is allowed. While this is not poor enough to reject by itself 
the hypothesis of no elongation, the fit allowing separate dependence describes the data better 
with a confidence level of 38%. When Eq. (132]) is used the fit not allowing separate dependence 
is rejected by a confidence level of 10" 7 , while allowing it results in the acceptable confidence 
level of 2%. For the fits of both Eq. (13T|) and Eq. (132]) the differences in x 2 between the fit 
allowing separate dependence and the fit not allowing it is huge, 296 and 464, respectively, 
while a difference of only 2 is expected if the hypothesis of no separate dependence is correct. 
This provides extremely strong evidence against the hypothesis of identical dependence of R2 
on the different components of Q. 

The fits using Eq. fl32l) are compared to the data in Fig.[TT] for small values of Q. As is 
apparent in Fig. [11] the shape of the dependence of R2 on the different components of Q is 
different at small values of these components and can not be described by equal values of R^, 

-Rsidc an d ?"out- 

We have repeated the fits varying the upper limit of Q to as low as 1 GeV. While some 
parameters show a considerable dependence on the Q-range of the fit, R S i de /RL and r out /R^ 
are found to be stable, varying by less than a standard deviation. The variation in the other 
parameters is thought to arise from correlations among the parameter values and the inability 
to determine well the values of the long-range parameters e^, e si d e and e out when Q is confined to 
small values. The value of the elongation, R sidc /R^, is far from unity and agrees well with the 
value in the fit of Eq. fj27]) imposing Eq. (1141) (Table [5]) and with that found using the Edgeworth 
parametrization. The value of r out /i?L is also far from unity, but is greater than unity, in 
contrast to Rside/Rh- Even more striking is the inequality of R S id c and r out : r out /-R S idc — 2. 
This indicates the invalidity of the azimuthal symmetry assumed in the elongation analyses. 

We also investigate elongation in the parametrization of Eq. (1111) by replacing AtQ 2 by 
ATiong^L + A r sidc<3sidc + ^ r out?ouf Fits of the resulting parametrization are performed both 
with Ariong, Ar S idc and Ar out as free parameters and with them constrained to be equal. In 
these fits only one m t bin is used, < m t < 4 GeV. The results are shown in Table [7] The 
conclusions are the same as for the fits of Eq. (130]) . The value of Ar si de/Ar L agrees with that of 
-Rside/ Rl, and Ar out is greater than unity. Constraining the fit by Atl = Ar sidc = Ar out results 
in an unacceptable \ 2 ■ 

We note that in all of the fits of ad hoc r-model parametrizations the long-range correlation 
parameters, e S id e and e out , deviate substantially and significantly from zero. This could indicate 
that not all non-BEC correlations are removed from the reference sample, which in turn could 
influence the amount of elongation found by the fits. To investigate this we have performed fits 
of the region Q > 1.8 GeV fixing the BEC contribution to zero (A = 0). The fits for Q < 4 GeV 
were then performed with the ej fixed to the values obtained in the A = fits. Of course, the 
confidence levels were lower than those of the fits where the e« could vary, but the elongation 
varied only slightly. For the fits of Table [6]^^/^ = 0.58±0.02 and R out /R^ = 0.984±0.002 
with a CL of 8%, while R sidc /R L = 0.64±0.02 and R out /R L = 1.25±0.03 with a CL of 0.09%. We 
also repeated the fits with the 6j fixed to zero. The confidence levels were much worse, but the 
elongation parameters proved quite robust: R S ide/ Rle = 0.69±0.02 and Rqui/Rle = 0.99±0.04 
with a CL of 10" 11 , and R sidc /RL = 0.72 ± 0.02 and R nt/Rh = 1-24 ± 0.03 with a CL of 10~ 13 . 
Further, note in Fig. [TT] that the data themselves show different dependencies on Ql, Qside, Qont- 
We conclude that the elongation is real, i.e., not an artifact of the parametrizations previously 
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Table 6: Results for two-jet events of fits of Eq. and Eq. (}3"Tj) or (|3"2]1 for Q < 4 GeV with 
all radii as free parameters and with them constrained to be equal. Uncertainties are statistical 
only. 



Eq.flHD 


A 


0.51 ±0.03 


0.49 ±0.03 




a 


0.46 ±0.01 


0.46 ±0.01 




Rle (fm) 


0.84 ± 0.04 


0.71 ±0.04 




Rside/ Rle 


0.60 ±0.02 


1 




Rout/ Rle 


0.986 ±0.003 


1 




e LE (GeV" 1 ) 


0.001 ±0.001 


0.000 ±0.001 




e S ide (GeV" 1 ) 


-0.069 ± 0.003 


-0.064 ±0.003 




e ut (GeV" 1 ) 


-0.032 ±0.002 


-0.035 ± 0.002 




7 


1.010 ±0.002 


1.012 ±0.002 




rVDoF 


14590/14538 


14886/14540 




CL 


38% 


2% 


Eq.fl32) 


A 


0.65 ±0.03 


0.57 ±0.03 




a 


0.41 ±0.01 


0.44 ±0.01 




Rl (fm) 


0.96 ±0.05 


0.82 ±0.04 




-Rside/ Rl 


0.62 ±0.02 


1 




r Q ut/RL 


1.23 ±0.03 


1 




€l (GeV- 1 ) 


0.004 ±0.001 


0.003 ±0.001 




e sidc (GeV" 1 ) 


-0.067 ±0.003 


-0.059 ± 0.003 




e out (GeV" 1 ) 


-0.022 ± 0.003 


-0.029 ± 0.002 




7 


1.000 ±0.002 


1.003 ±0.002 




^/DoF 


10966/10647 


11430/10649 




CL 


2% 


lO" 7 
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Figure 11: Projections of the Bose-Einstein correlation function R2(Ql, Qside, Qout) for two-jet 
events using the region below 80 MeV for the non-projected components. 
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Table 7: Results for two-jet events of fits of Eq. (TJTJ with ArQ 2 replaced by Ari on g<3L + 
Ar side (5g ide + Ar out go Ut and long-range parametrization 1 + e L Q L + e sidc <3 sidc + e out g out , where e L , 
e sido and e out are free parameters. In the first fit all three Ar parameters are free; in the second 
they are constrained to be equal. Uncertainties are statistical only. 



A 


n c7 _i_ n no 


n a o _i_ n no 

u.4y ± U.Uz 


a 


0.47 ±0.01 


0.50 ±0.01 


Ar L (fm) 


1.74 ±0.15 


1.32 ±0.10 


Ar sidc /Ar L 


0.61 ±0.03 


1 


Ar out /Ar L 


1.48 ±0.08 


1 


e L (GeV- 1 ) 


0.000 ± 0.001 


-0.001 ±0.001 


e S idc (GeV- 1 ) 


-0.064 ± 0.003 


-0.056 ±0.002 


eout (GeV- 1 ) 


-0.021 ±0.003 


-0.025 ±0.002 


7 


1.008 ±0.002 


1.010 ±0.002 


X 2 /DoF 


10977/10647 


11239/10649 

io- 5 


CL 


1% 



used, and that the ad hoc modified r-model provides a reasonable description of the elongation. 

6 Discussion and Conclusions 

The usual parametrizations of BEC of pion pairs in terms of their four-momentum difference Q, 
such as the Gaussian, Edgeworth, or symmetric Levy parametrizations, are found to be inca- 
pable of describing the Bose-Einstein correlation function R 2 , particularly the anti-correlation 
region, approximately 0.5 < Q < 1.5 GeV. This failure has not been previously so obvious, since 
previous analyses generally fit R 2 (Q) only up to 2 GeV or less and the anti-correlation region is 
then tacitly absorbed into the parametrization of long-range correlations. The existence of this 
anti-correlation region has the more general implication of ruling out the form usually proposed 
for BEC, R 2 = l + \f(Q)\ 2 . 

Those parametrizations of BEC are based on a static view of pion emission. No time 
dependence is assumed, which means that either the pion emission volume is unchanging during 
pion emission or that the parameters describing the volume, which result from fitting R 2 , are 
some sort of time average. 

The r-model assumes a very high degree of correlation between momentum space and space- 
time and introduces a time dependence explicitly. The Bose-Einstein correlation function R 2 
measures the real part of products of Fourier-transformed proper-time distributions. Not only 
the positive correlations at low values Q but also negative correlations at somewhat higher 
values of Q are predicted. It is as though the Bose-Einstein symmetry pulls identical pions 
below some critical value of Q closer together creating an excess of pairs at lower Q and leaving 
a deficit of pairs around this value of Q. 

In the r-model R 2 is then a function of one two-particle variable, the invariant four- 
momentum difference of the two particles Q, and of two single-particle variables, the values 
of a of the two particles, where a is a parameter in the correlation between momentum space 
and space-time, Eq. ([2]), The time-dependence is assumed to be given by an asymmetric Levy 
distribution, H(t), which imposes causality by allowing particle production only for r > 0. 
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The parameter a can be absorbed into an effective radius to obtain an expression, Eq. ffTB"]) . for 
R2{Q) which successfully describes BEC, including the anti-correlation region, for both two- 
and three-jet events. We note that a similar anti-correlation region has recently been observed 
by the CMS Collaboration [TU] in pp collisions at y/s = 0.9 and 7TeV and successfully fit by 
the r- model formula, Eq. ffT3~j) . 

For two-jet events a = l/m t and the introduction of an effective radius is unnecessary. The 
BEC correlation function is then a function not only of Q but also of the transverse masses of 
the pions. This description of R2(Q,m t i,m t2 ) also successfully describes the data. 

Nevertheless, the r-model description breaks down when confronted with data expressed 
in components of Q. The r-model predicts that the only two-particle variable entering R 2 
is Q, i.e., that there is a single radius parameter which is the same for each of the separate 
components of Q. This implies a spherical shape of the pion region of homogeneity, whereas 
previous analyses using static parametrizations have found a shape elongated along the event 
axis, as evidenced by -R S ide being smaller than in the LCMS. Assuming azimuthal symmetry 
about the event axis, i? S ide is the transverse radius of the region of homogeneity. 

Accordingly, the r-model equations for R 2 have been modified ad hoc to allow an elongation. 
Moreover, fits have been performed not only in the LCMS, but also in the rest frame of the pair. 
Fits not allowing elongation have much worse confidence levels than fits allowing elongation. 
The elongation found agrees with that found in the previous analyses. However, in the rest 
frame of the pair the radius parameter in the out direction is found to be approximately twice 
that in the side direction. Thus the assumption of azimuthal symmetry about the event axis is 
found not to hold. 

We note that an elongation is expected [80J for a hadronizing string in the Lund model, 
where the transverse and longitudinal components of a particle's momentum arise from different 
mechanisms. Perhaps the absence of azimuthal symmetry arises from gluon radiation, which is 
represented by kinks in the Lund string. Therefore a possible improvement of the r-model could 
be to allow the longitudinal and transverse components to have different degrees of correlation 
between a particle's momentum and its space-time production point in Eq. (|8~|) . 

The r-model, as we have used it, incorporates a Levy distribution. A Levy distribution 
arises naturally from a fractal, or from a random walk or anomalous diffusion [71], and the 
parton shower of the leading log approximation of QCD is a fractal [SBHTO] . In this case, the 
Levy index of stability is related to the strong coupling constant, a s , by [511152"] 

a s = —a . (33) 

Assuming (generalized) local parton hadron duality [83H85] . one can expect that the distribution 
of hadrons retains the features of the gluon distribution. For the value of a found in the fit of 
Eq. (HI]) for two-jet events (Table E} we find a s = 0.46 ±0.01 

-ao2- This is a reasonable value for 
a scale of the order of 1 GeV, which is the scale at which the production of hadrons is thought 
to take place. For comparison, from r decay, a s (m T « 1.8 GeV) = 0.34 ± 0.03 [SB], and from 
the average value of a s (M z ), a s (l GeV) = 0.50±g;g| [UJ. 

It is of particular interest to point out the m t dependence of the "width" of the source. 
In Eq. (fiT|) the parameter associated with the width is At. Note that it enters Eq. (TIT!) as 
AtQ 2 /m^. In a Gaussian parametrization the radius R enters the parametrization as R 2 Q 2 . 
Our observation that At is independent of m t thus corresponds to R oc l/y/W t and can be 
interpreted as a confirmation of the observation [20TI22] of such a dependence of the Gaussian 
radii in 2- and 3-dimensional analyses of Z decays. The lack of dependence of all the parameters 
of Eq. (Hip on the transverse mass is in accordance with the r-model. 
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Given these successes, the BEC fit results of the r-model have been used to reconstruct the 
pion emission function of two-jet events. Note that all charged tracks are treated as pions. For 
BEC the main influence of this is a lowering of the value of the parameter A. However, the 
influence of this as well as the influence of resonances on the rapidity and p t distributions may 
be more profound. Hence, the emission function is subject to some quantitative uncertainty. 
Nevertheless, the general features of the emission function should remain valid: Particle pro- 
duction begins immediately after collision, increases rapidly and then decreases slowly, occuring 
predominantly close to the light cone. In the transverse plane a ring-like structure expands out- 
wards, which is similar to the picture in hadron-hadron interactions but unlike that of heavy-ion 
collisions [5]. Despite this similarity the physical process is much different. Reflecting the non- 
thermal nature of e + e~ annihilation, the proper-time distribution and the space-time structure 
are reconstructed here without any reference to a temperature, in contrast to the results of 
earlier hadron-hadron and heavy-ion collisions. 
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Appendix 

In Tables [8] and M we present the data used in the fits of Section [3j 



Table 8: Data for fits to B,2(Q) for two-jet events used 
in Figs. El HI and El as well as Tables (H El and |3j Uncer- 
tainties are statistical only. 



Q (GeV) 


(mti) (GeV) 


(m t2 ) (GeV) 


R2 


0.030 


0.265 


0.245 


1.362±0.051 


0.064 


0.283 


0.238 


1.424±0.021 


0.102 


0.303 


0.232 


1.362±0.013 


0.141 


0.320 


0.225 


1.270±0.009 


0.181 


0.337 


0.222 


1.210±0.007 


0.221 


0.354 


0.219 


1.154±0.006 


0.260 


0.370 


0.218 


1.118±0.005 


0.300 


0.387 


0.220 


1.091±0.005 


0.340 


0.402 


0.222 


1.062±0.005 


0.380 


0.418 


0.225 


1.028±0.004 


0.420 


0.432 


0.229 


1.023±0.004 


0.460 


0.445 


0.233 


1.005±0.004 


0.500 


0.457 


0.238 


0.991±0.004 


0.540 


0.472 


0.243 


0.979±0.004 


0.580 


0.485 


0.247 


0.968±0.004 


0.620 


0.496 


0.252 


0.965±0.004 


0.660 


0.508 


0.256 


0.962±0.004 



Continued on next page 
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Table 8 — continued from previous page 



Q (GeV) 


(m tl ) (GeV) 


(m t2 ) (GeV) 


R2 


0.700 


0.519 


0.261 


0.959±0.004 


0.740 


0.531 


0.265 


0.959±0.004 


0.780 


0.540 


0.268 


0.959±0.004 


0.820 


0.549 


0.271 


0.958±0.004 


0.860 


0.559 


0.275 


0.957±0.005 


0.900 


0.565 


0.278 


0.963±0.005 


0.940 


0.572 


0.281 


0.958±0.005 


0.980 


0.579 


0.284 


0.965±0.005 


1.020 


0.582 


0.285 


0.962±0.005 


1.060 


0.589 


0.287 


0.967±0.005 


1.100 


0.594 


0.290 


0.965±0.005 


1.140 


0.599 


0.292 


0.962±0.005 


1.180 


0.602 


0.293 


0.969±0.005 


1.220 


0.605 


0.295 


0.962±0.005 


1.260 


0.611 


0.296 


0.988±0.006 


1.300 


0.612 


0.298 


0.982±0.006 


1.340 


0.617 


0.298 


0.985±0.006 


1.380 


0.616 


0.300 


0.979±0.006 


1.420 


0.619 


0.301 


0.979±0.006 


1.460 


0.621 


0.301 


0.992±0.006 


1.500 


0.625 


0.302 


0.975±0.006 


1.540 


0.622 


0.303 


0.988±0.007 


1.580 


0.623 


0.304 


0.990±0.007 


1.620 


0.629 


0.305 


0.983±0.007 


1.660 


0.627 


0.304 


0.991±0.007 


1.700 


0.630 


0.306 


0.995±0.007 


1.740 


0.631 


0.307 


0.984±0.007 


1.780 


0.629 


0.306 


0.990±0.007 


1.820 


0.636 


0.308 


0.982±0.007 


1.860 


0.632 


0.307 


0.993±0.007 


1.900 


0.635 


0.309 


0.979±0.007 


1.940 


0.634 


0.308 


0.986±0.008 


1.980 


0.635 


0.309 


0.991±0.008 


2.020 


0.634 


0.311 


0.995±0.008 


2.060 


0.631 


0.309 


0.998±0.008 


2.100 


0.636 


0.310 


0.993±0.008 


2.140 


0.634 


0.310 


0.990±0.008 


2.180 


0.637 


0.311 


0.998±0.008 


2.220 


0.636 


0.312 


0.985±0.008 


2.260 


0.638 


0.312 


1.000±0.009 


2.300 


0.635 


0.311 


0.975±0.008 


2.340 


0.635 


0.312 


0.998±0.009 


2.380 


0.636 


0.313 


0.998±0.009 


2.420 


0.634 


0.312 


0.978±0.009 


Continued on next page 
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Table 8 — continued from previous page 



Q (GeV) 


(m tl ) (GeV) 


(m t2 ) (GeV) 


R2 


2.460 


0.636 


0.311 


0.986±0.009 


2.500 


0.637 


0.314 


0.985±0.009 


2.540 


0.635 


0.312 


1.000±0.009 


2.580 


0.636 


0.312 


0.989±0.009 


2.620 


0.635 


0.313 


0.988±0.009 


2.660 


0.633 


0.312 


0.999±0.010 


2.700 


0.635 


0.311 


0.981±0.010 


2.740 


0.635 


0.313 


0.983±0.010 


2.780 


0.633 


0.311 


1.001±0.010 


2.820 


0.632 


0.312 


1.003±0.010 


2.860 


0.633 


0.312 


1.006±0.010 


2.900 


0.633 


0.313 


0.981±0.010 


2.940 


0.631 


0.310 


1.010±0.011 


2.980 


0.632 


0.311 


1.018±0.011 


3.020 


0.629 


0.310 


0.993±0.011 


3.060 


0.634 


0.312 


0.990±0.011 


3.100 


0.631 


0.310 


1.005±0.011 


3.140 


0.635 


0.310 


0.983±0.011 


3.180 


0.634 


0.313 


0.991±0.011 


3.220 


0.635 


0.311 


0.998±0.011 


3.260 


0.626 


0.309 


1.001±0.011 


3.300 


0.633 


0.312 


0.994±0.011 


3.340 


0.634 


0.312 


0.981±0.011 


3.380 


0.633 


0.310 


0.980±0.011 


3.420 


0.632 


0.311 


0.983±0.012 


3.460 


0.630 


0.310 


0.982±0.012 


3.500 


0.632 


0.310 


0.989±0.012 


3.540 


0.636 


0.311 


0.990±0.012 


3.580 


0.633 


0.311 


0.996±0.012 


3.620 


0.629 


0.310 


1.001±0.012 


3.660 


0.633 


0.312 


0.999±0.012 


3.700 


0.633 


0.311 


0.992±0.012 


3.740 


0.635 


0.310 


0.988±0.012 


3.780 


0.635 


0.313 


0.999±0.013 


3.820 


0.636 


0.312 


0.986±0.013 


3.860 


0.637 


0.309 


0.997±0.013 


3.900 


0.631 


0.309 


0.979±0.013 


3.940 


0.633 


0.309 


0.992±0.013 


3.980 


0.636 


0.313 


1.009±0.013 
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Table 9: Data for fits to R 2 {Q) for three-jet events used in Figs. [3] and [5] as well as Tables [T] 
and El Uncertainties are statistical only. 



Q(GeV) R 2 Q(GeV) R 2 Q (GeV) R 2 
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